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Abstract 

The minimal Standard Model extension with the Weinberg operator does accommodate the 
observed neutrino masses and mixing, but predicts a neutrinoless double beta (Oi//3/3) decay rate 
proportional to the effective electron neutrino mass, which can be then arbitrarily small within 
present experimental limits. However, in general Ov/SlB decay can have an independent origin 
and be near its present experimental bound; whereas neutrino masses are generated radiatively, 
contributing negligibly to Oiyf3f3 decay. We provide a realization of this scenario in a simple, well 
defined and testable model, with potential LHC effects and calculable neutrino masses, whose two- 
loop expression we derive exactly. We also discuss the connection of this model to others that 
have appeared in the literature, and remark on the significant differences that result from various 
choices of quantum number assignments and symmetry assumptions. In this type of models lepton 
flavor violating rates arc also preferred to be relatively large, at the reach of foreseen experiments. 
Interestingly enough, in our model this stands for a large third mixing angle, sin 

6*13 Z 0.008, when 

/i — > eee is required to lie below its present experimental limit. 



1 



I. INTRODUCTION 



Neutrino oscillations are the only new physics (NP) beyond the minimal Standard Model 
(SM) observed to date [1] (for recent reviews see for instance [2-4]). They can be fully 
explained introducing rather small neutrino masses m,, and the corresponding (unitary) 
charged current mixing matrix U [5, 6]. The observed pattern of neutrino masses can be 
implemented, however, in two quite different ways depending on the Dirac or Majorana 
character of the neutrinos. In the first case the SM is extended by adding three SM singlets, 
to provide Dirac masses to the SM neutrinos, Ui, through small Yukawa couplings. 
Alternatively, we can consider extending the SM by adding the only invariant dimension 5 
(Weinberg) operator that can be written using the SM field content [7] ^, O5 = C.L<p(j)^ ii- In 
this case the SM neutrinos acquire Majorana masses after electroweak symmetry breaking, 
(0°) = and are inversely proportional to A, the NP scale associated with O^. The small 
neutrino masses then require that the coefficients of this operator be small, either due to a 
very large NP scale or suppressed dimensionless couplings. 

Both alternatives (that light neutrinos are Dirac or Majorana), are viable and indistin- 
guishable if A is very large, except for the possible observation of lepton number violation 
(LNV) in neutrinoless double beta (OuPP) decay [8] ^ (for a recent review see [12]). Indeed, 
the SM extension with Dirac neutrino masses preserves lepton number (LN), and hence OuPP 
decay is forbidden. Whereas Majorana masses carry LN equal two, as does the Weinberg 
operator, and this leads to a non-zero Oz//3/3 width. Thus, the observation of OuPP decay 
would strongly favor Majorana neutrino masses [13]; hence the prime relevance of this type 
of experiments (see [14, 15] for recent reviews). 

All this, however, relies on the assumption that both such minimal SM extensions describe 
the dominating NP effects up to very high scales. We will argue, however, that LNV and 
neutrino masses may be due to NP near the electroweak scale, in which case a much richer set 
of possibilities can be realized. The main point we wish to emphasize is that although both 
OuPP decay and Majorana neutrino masses do violate LN, they need not be both directly 
related to the Weinberg operator O^, unlike the above minimal SM extension consisting 

^ (p and £l are the SM Higgs and lepton doublets and (j> — iT2(p* , £l — *'''2^i- 

^ We shall not discuss the implications of requiring enough leptogenesis to account for the observed baryo- 
genesis [9] (for recent reviews see [10, 11]). 
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only of this dimension 5 operator. For instance, the leading LNV effects in Ou/S/S decay 
could be mediated by an operator involving only right-handed (RH) electrons, such as, for 
example, Og = (e^e^) (^(p'^D^cfj (^(p'^D^cfj which is the lowest-order LNV operator with two 
RH leptons and invariant under the SM gauge transformations. In this case O5 (with the 
associated Majorana masses) is generated radiatively by Og at two loops, where the charged 
leptons suffer a chirality change through mass insertions and are transformed into neutrinos 
through the exchange of a boson. This results in a further suppression by two charged 
lepton masses divided by two powers of the NP mass scale, which we do assume to be close to 
the electroweak scale. These radiatively-generated Majorana masses will produce the usual 
contribution to Oz//3/3 decay, which, however, will be negligible compared to the Og one. 

In a companion paper [16] we classify the different ways of generating Oz//3/3 decay and 
light neutrino masses by the addition of higher order effective operators. This has been 
studied in the literature [17-20], but mostly for operators involving fermions and scalars; 
we will concentrate instead on operators involving gauge bosons but not quarks (thus, for 
example, excluding ab initio models with heavy leptoquarks from our analysis). Here we 
will instead provide a realistic testable model realizing the above scenario, where (i) LN is 
broken at the electroweak scale; (ii) Oz//3/3 decay into two RH electrons has a rate of the 
order of its experimental limit, through the tree-level exchange of new scalars; and (Hi) it 
contains finite, and therefore calculable, neutrino masses. Despite a relatively small number 
of parameters this model can also accommodate the observed pattern of neutrino masses 
and mixings, which are generated at two-loop order (and whose contribution to OuPP decay 
is in this case negligible). This model is related to others that have been discussed in the 
literature [21, 22], the differences are crucial, and essential in maintaining the 3 features just 
mentioned; we discuss these points in detail below. 

There are many SM extensions where Ou/SP decay receive new contributions besides those 
from light Majorana neutrino masses (see [23] for a general overview). The simplest scenario 
assumes the presence of heavy Majorana neutrinos whose exchange generates new contribu- 
tions to Oz//3/3 decay, similar to the ones generated by the light neutrinos. (See for recent 
work [24-27].) Other extensions with many more new particles have also been studied; such 
as left-right (LR) models [28-30] and supersymmetric extensions (see, for instance, for a re- 
view [1], and references there in). In such models there are several contributions to the Oz//3/3 
decay amplitude, any of which can dominate depending on the region of parameter space 
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being considered. We are interested, however, in identifying the minimal SM extensions 
leading to the largest possible contributions to Oz//3/3 decay and containing no independent 
neutrino masses. This means that, as argued above, light neutrino masses only result from 
the unavoidable contributions mediated by the LNV operators generating Ou/SP decay; so 
that the effective Lagrangian approach provides the proper language for classifying such 
scenarios [16]. Here, we are only concerned with giving a specific example of the case where 
OuPP and neutrino masses are generated by the same underlying physics but at different 
orders in perturbation theory, with the former appearing at tree level, while the latter only 
at two loops ^. Simple models with finite neutrino masses at two loops have been often 
discussed in the literature [32, 33], although not necessarily related to NP inducing Ou/SP 
decay as it is the case here. 

It is worth emphasizing again that the model we study is one particularly simple example 
of a class of models realizing the above scenario, and that all such models share many 
phenomenological implications. In our case the model contains a discrete Z2 symmetry, 
which is spontaneously broken. This model has the virtue of providing a direct analysis 
of the symmetries and scales, however, it also has a domain- wall problem [34, 35]. One 
way of dealing with this problem is to allow the discrete symmetry to remain unbroken, 
but at the price of generating both neutrino masses and the OuPP amplitude at a higher 
loop order; a possibility we will not pursue. We instead discuss related but somewhat more 
involved models that avoid the domain wall problem while retaining the same low-energy 
phenomenology. We will restrict ourselves to the most relevant region of parameter space 
where the OuPP decay rate lies within the reach of the next round of experiments. With this 
assumption, together with the constraints from lepton flavor violation (LFV) processes such 
as /i — )■ eee, and with the requirement of perturbative unitarity (or if preferred, naturalness of 
perturbation theory), the model predicts that the neutrino masses obey a normal hierarchy, 
that the lightest neutrino mass lies in the range 0.002 eV ^ mi ^ 0.007 eV, and that the 
third mixing angle in the Pontecorvo-Maki-Nakagawa-Sakata mixing matrix U [5, 6] obeys 
sin^6'i3 ^ 0.008. Value which lies well within the sensitivity of ongoing neutrino oscillation 

^ LNV effective operators including quarfcs also generate neutrino masses but in general at higher loop 
order, in fact too high in some cases to explain the observed spectrum [31]. Although such operators are 
not considered in our analysis, we will comment on them further when discussing our general set up in 
detail [16]. 
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experiments. Besides, the new (charged) scalar masses can be within the LHC reach, and 
various LFV processes are predicted to have rates that will be probed when present precision 
is improved by the next generation of experiments. 

In next section we present our model; few details on the scalar spectrum and couplings 
are also worked out, and its relation to other models is briefly discussed. The reader mainly 
interested in the phenomenological implications of the model can go directly to Section 
III, where we evaluate the rate for Oz//3/3 decay. The requirement that this process can 
be observed in the next round of experiments, together with perturbativity, translates into 
upper bounds on the masses of the extra scalars. These, in turn, result in limits on their 
couplings if present bounds on LFV processes are to be fulfilled, as we show in Section IV. We 
calculate the neutrino masses in Section V, where we show that the model can accommodate 
the observed pattern of neutrino masses and mixings, though it predicts a somewhat small 
electron neutrino effective mass mee- Finally, the prospects for the discovery of the extra 
scalars at LHC are considered in Section VI. In connection with this it must be noted 
that models with extra scalars may allow for a would-be Higgs boson more elusive to LHC 
searches, for it can have further decay channels open (see for instance [36-38]). Conclusions 
are drawn in last section. In three Appendices we collect some technical details. 

II. A MODEL WITH LEPTON NUMBER SOFTLY BROKEN 

This model only extends the SM Higgs sector, in order to allow for scalar couplings to 
lepton bilinears with non-zero LN. More precisely, as we look for separating the origin of 
neutrino masses from the mechanism mediating Ou/SP decay, we only introduce new scalar 
couplings to RH charged leptons, ensuring they are the only final state produced in tree-level 
Oi//3/3 decay. A simple way to achieve this is to enlarge the Higgs sector including, besides 
the SM scalar doublet of hypercharge 1/2, a complex scalar singlet k, of hypercharge 2, 
a real (neutral) scalar singlet a, and an electroweak triplet x of hypercharge L We also 
impose a Z2 symmetry, under which all SM particles and k, are even while both cr and x are 
odd, which forbids their coupling to lepton bilinears. Thus, the Yukawa Lagrangian reads 

Cy = iLYeCR (f) + e'j^geRK + h.c. , (1) 
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where we can assume Ye to be a 3 x 3 diagonal matrix with positive eigenvalues, and g a 
complex symmetric matrix with three of its phases unphysical. 

The most general Higgs potential consistent with the symmetries is ^ 

+ ^<r^' + ^M' + Ax (Tr {x^x]Y + A^Tr { {x\f] 
+A^.|0| V + X^M?^? + A^xl^l'Tr {xh] + A;^0tx^x0 
+A<^«(tVP + X^^\k\^T:i [x^x] + A^^^or^Tr [x^x] 

+ [/i««:Tr|(x^)'} + A6a0t;^0 + h.c.] , (2) 

where by rephasing the fields we can always choose /x^ and Ag real of either sign. For 
convenience, we will take /i^ positive and \ negative. In terms of charge eigenstates, the 
fields X a-iid 4' ^"^^ written 




X= ' _ , 0= J • (3) 



The singlet a is introduced to preserve the discrete symmetry that forbids the scalar triplet 
from coupling to leptons. This symmetry is broken spontaneously by the VEVs (cr) and 

For the phenomenological discussion it is important to note that in the limit of vanishing 
lepton number can be defined to only act on this is enough to protect neutrinos from 
getting a Majorana mass. While the lepton number, which would forbid Oz//3/3 decay 
(into two e^j), is explicitly broken in the scalar potential due to the presence of the and 
Ae terms. Therefore, the Oz//3/3 decay amplitude will be proportional to all three couplings: 
Qab in Eq. (1) and and Ae in Eq.(2); whereas neutrino masses will also depend on Yg- 

This is one of a class of models with the same low energy physics. In order to understand 
the common LNV features it is convenient to consider the above theory, but with a com- 
plex, and the modifications needed to insure a real Lagrangian. In this case LN is exactly 
conserved; if leptons are assigned LN —1, k carries LN equal to 2, while cr^ and x both 
carry LN equal to 1, so that type II see-saw Yukawa couplings iiX^L are forbidden. The 



^ Terms like Tr {x^} Tr {x^^}, Tr {x^X^^} or 4>^XX^4' can easily be related to the terms already included by 
using the fact that for two arbitrary 2x2 traceless matrices A = AiTi and B = BiTi, {A, B} = Tr {A_B} 1. 
Other possible terms are forbidden by the discrete symmetry. 
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vertex ctHlX^l, however, is allowed and generated at two loops, providing light neutrinos a 
finite mass (which we calculate exactly) after spontaneous symmetry breaking proportional 
to (o"^x)- Obviously, this model has a Majoron [39] because LN, which is a global abelian 
symmetry, is spontaneously broken. One can then consider the region of parameter space 
where such interesting models are phenomenologically viable ^, or promote the symmetry 
to a local one by gauging, for instance, baryon minus lepton number (B— L). Although, this 
is an interesting possibility, too, with renewed experimental interest (see, for instance, for 
a review [45]), it requires adding RH neutrinos that provide a new potential source of light 
neutrino masses; this lies outside the goal of the present investigation that is centered on 
theories where neutrino masses are generated by the LNV effective operator accounting for 
Ou/Sp decay. 

For this reason we have assumed cr to be real, breaking LN explicitly but leaving as a 
remnant an exact Z2 symmetry. The general features of the theory remain, though the 
Majoron does not appear and LN is reduced to the above Z2 symmetry ^. The spontaneous 
breaking of this discrete symmetry will generate domain walls. This poses no problem 
provided (cr) is sufficiently large (with Ag correspondingly small while keeping their product 
fixed) guaranteeing the formation of such defects before the infiationary epoch. We have 
assumed instead that (cr) is not too large because we prefer avoid requiring Ae unnaturally 
small; we may then avoid a domain-wall problem by adding soft-breaking terms, such as 
to the potential (such a modification would not require the introduction of tree-level Yukawa 
interactions violating the discrete Z2 symmetry). We do not pursue this discussion because 
within this class of models there is an even simpler one, with the same neutrino physics at 
low energy, and none of these potential drawbacks. It is the same model presented above, 
but with a replaced by its vacuum expectation value, a — > (cr); up to coupling constant 
redefinition this yields the same potential as Eq. (2) without the terms containing cr, except 

^ In this model the Majoron will be mainly a singlet [39] . In which case its couplings to ordinary matter are 
small and then little constrained, while its coupling to the Higgs boson is essentially free. Thus, singlet 
Majorons can result in invisible Higgs decays (see [40] for a recent example). Moreover, they can have 
interesting implications in astrophysics and cosmology because they can substantially affect the cooling 
of supernovas [41] and the neutrino relic abundance, due to the possibility of neutrino decaying [42] or 
annihilating [43] into Majorons. They can be even used as a dark matter candidate (see for instance [44]) 
if massive (pseudo-Majorons). 

^ One can also assign an odd Z2 parity to the leptons, in which case this discrete symmetry equals (— 1)^^. 
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for the last term that becomes fi^(l)^x4'j with 

f^x = ^6 {(t) , (4) 

of the order of a TeV. Obviously, the resulting renormalizable Lagrangian has the same 
quantum behaviour than ours; in particular, neutrino masses are finite and generated at two 
loop order and can be obtained from our results by eliminating Xq using (4). We will prefer 
to discuss the model including the real scalar field cr because the perturbative and symmetry 
analyses appear to be more transparent to us. In the minimal model with the SM addition 
of only K and and the scalar potential terms relating their LN charges ^hkTi | 
f^x^^X^P the LN of the x field is not well defined since the first term requires it to be 1 while 
the second (though never 2); as a consequence the effective vertex Y^i^x^L, is finite and 
generated radiatively Despite the parallels of this discussion with models implementing a 
type II see-saw mechanism (which contain a tree- level IlX^l coupling) there is an important 
difference, namely, that in such theories both the expression for the neutrino masses and the 
Oi//3/3 decay amplitude are linear in (%); in contrast the corresponding results in our model 
are proportional to (x)^ (see Eqs. (25) and (33) below). 

The extension of the SM model we consider is related to the one presented in Refs. [21, 22] 
as far as particle content is concerned (we differ by adding the singlet cr). However, the 
symmetries and quantum number assignments are different, which proves a crucial difference. 
Had we included the hard term K.(f)^x^'P &s was done in the above references, x should have 
been assigned LN = 2; this would necessitate also including the tree-level coupling iiX^L 
that would lead to the usual type II see-saw scenario. In particular it would have been 
inconsistent to assign LN zero to x oi' fo arbitrarily exclude this Yukawa coupling from the 
Lagrangian, for its coefficient would receive divergent radiative corrections; in consequence 
the neutrino masses are not calculable. If one requires x fo have LN equal to 0, as was done 
in these publications, the quartic coupling must be absent; but then LN remains 

unbroken and the light neutrino masses must vanish to all orders which is again inconsistent 
with the results presented there. 

Note that in the theory containing an extra real scalar a, one could also reason differently to justify this 
result, fndeed, we could assign cr LN equal to 0; then the quartic term a<j>^x4' fixes the x LN also equal to 
0, while the trilinear term KTr{(x^^)^} breaks LN softly. In any case the discrete Z2 symmetry guarantees 
that the neutrino masses stay finite. 



8 



A. The scalar spectrum 



The requirement on the scalar potential of being bounded from below is fulfilled restricting 
the quartic couplings in Eq. (2) adequately; these conditions include ® 

•^cr; > , A^CT > •~2a/Ao-A^ . (5) 

We have also checked that there is a non-trivial minimum on which the scalar neutral 

components acquire non-zero expectation values: = {(fP) > 0, = > 0, i;^ = {a) > 
0, with 

/ = (/> + + ^'^^) ' = (X°> + -^{XR + ^Xi) , cT={a) + an . (6) 

In Appendix A we comment on the experimental limit on the scalar triplet VEV f^, which 
we will find to be of the order of few GeV; to be conservative we will assume f ^ < 5 GeV, 
this satisfies f ^ ^ f <^ 174GeV where l/(t'| + 2f^) = g'^/2m'^ = 2\f2Gp (derived from Eq. 
(A2) and the experimental limit on the p parameter, \pq — 1| ^ 1 [1]). In this approximation 
the minimization conditions can be easily solved: 

~ I_ q)^ ~ I_ -ll^ ~ Z_ (Y) 

" ml + vlX^^ + vlK^ ' ^ 4A.A^ - A|, ' " 4A.A^ - A^^ ' ^ ^ 

Notice that the phase choice Ag < is consistent with being real and positive. We set 
the a and (f) mass squared terms in the Higgs potential negative to favour the development 
of such a minimum. Though we choose the x niass squared term positive, this field also 
acquires a small VEV induced by the doublet and singlet VEVs, similar to the case observed 
in see-saw of type II models [46-48]. 

The scalar masses can be obtained by substituting Eq. (6) in the potential. Using the 
exact minimization conditions to eliminate and mg. in favour of the VEVs, the mass 



Notice that the term |(/)| cr is always positive 
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terms for the charged scalars can be written 



Ml 



Ml 



K 



++ 



++ 




, with 



2 \/ 



2vl 



2v^ 



-V2V^Vy^ \ 



(8) 
(9) 
(10) 



Analogously for the neutral sector 



Cm — 



{h Xi)Mf 




{(t>R Xr (7r)M^ 



Xr 



with 



VaXe 



(11) 



(12) 



M 



R 



-2v^v^ 

^ 4v|A</, 2v^ {vcrXe + Vy^X^^) 2^v^ [v^X^a + v^Ae) ^ 

2v^ {v^Xq + Vy^X^,^) 4v^(A^ + A^) - vIv^Xq/v-^ \f2 (u^Ae + 2vaVyXa^ (-13) 

2\f2v^ {v^X^^ + v^Xe) \f2 (w^Ag + 2v^v^X„^ Sv^A^, - 2v\v^XqIv„ ^ 

All eigenvalues of these mass matrices must be positive (except for the would-be Goldstone 
bosons providing the longitudinal vector boson degrees of freedom) in order to guarantee 
that the solution to the minimization conditions corresponds to a local minimum. This sets 
further constraints on the model parameters, which can be satisfied rather easily, especially 



in the limit ^ 



Thus, we are left with two massive doubly-charged scalars k\^2i 

K,2 — — sin^i)^"'""'" -I- cos^dX''"'' , 



with 



sin 29d = 2 sin 9d cos 9d 



ml, - ml^ 



and only one massive, mainly triplet, singly-charged scalar u, 

a;+ = -sin^50+ -Fcos^sX"^ ' "^^^^ tan 9s = "^^^ 



(14) 



(15) 



(16) 
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Similarly, there is a neutral scalar A with imaginary components, 

2v 

A = -sin6i(j)j + cosdixi > with tandj = — - . (17) 

There are also three neutral scalars along the real components (p^, xRi (^R- We will denote 
these mass eigenfields by h (mainly doublet), H (mainly triplet) and s (mainly singlet). 
They are obtained rotating the current fields, what introduces other three mixing angles. 
Notice that in the limit ^ f^^^o-, we have <^ v^^ri and mf52,tj,A ~ with all mixings 
small. 

B. Some scalar couplings of phenomenological interest 

Once the quadratic terms of the Lagrangian are diagonalized we can read the interactions 
for the mass eigenfields. In the following we will need the scalar coupling to RH electrons 

dRa 9ab em (cos6'i5ft;i - sin6'£,K2) + h.c. , (18) 

and the corresponding doubly-charged scalar couplings to gauge bosons 

gh^^W'W^-X^'^ + h.c. ^ g^W-W"- {smODKi + cos^d/«2) + h.c. , (19) 

as well as their trilinear couplings 

- {x-y - 2fi^K++x-x'^ + AeCT (0^)' + h-c , (20) 

which can be also expressed in terms of the mass eigenfields using Eqs. (14-16), and the 
corresponding VEVs in Eq. (6). Finally, the Yukawa coupling changing charge and chirality 
writes 

W Ye en 0+ + h.c. = VZYeCR (cos esG+ - sin ^50;+) + h.c. , (21) 

where = cos 6*50+ + sin 9sX^ is the would-be Goldstone boson providing the third com- 
ponent to the W. 

III. NEUTRINOLESS DOUBLE BETA DECAY 

Both doubly-charged scalars Ki,2 have components along the singlet and the triplet 
Therefore, they (respectively) contain couplings to RH electrons and to W^s, gen- 
erating an effective vertex eeWW that mediates 0^(3(3 decay. In this section we calculate 
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Figure 1: Dominant tree- level diagram contributing to the effective neutrinoless double beta decay 
operator. 

this contribution to Oz//3/3 decay and obtain the exphcit constraints on the model param- 
eters derived from the assumption that Oz//3/3 decay will be observed in the next round of 
experiments. 

Assuming that m^^^ 3> f 0,0- and integrating out the heavy k and x modes we find, after 
a straightforward calculation, that the effective Lagrangian contains the term 

= ^i^flS/^ (e]^ gl, e%,) U^D'^4>) U'D^^) + h.c. , (22) 

as announced in the Introduction, and discussed in the companion paper [16]. One can 
better understand the origin of this LNV interaction by considering the contribution of the 
dominant diagram in Fig. 1, where the different couplings and VEVs involved are displayed 
explicitly. The corresponding eeWW vertex at low energy (g^ <^ m^^ ^) can be written as ^ 

^-lyiy = - ^2\J ^a git e%, + h.c. , (23) 

where we have summed up all possible mass insertions in the internal propagator, and used 
Eq. (15). This expression coincides with Eq. (22) when the scalar doublets develop a VEV 
in the limit of large m^^^ (:» "^^0,0-)) ^-C-, ~ — Agfo-fl/m^ and ^ m^^, ^ m^a- 

Let us particularize to the case a = b = e, relevant for Oz//3/3, and further integrate the 

^ It must be emphasized that the Lagrangian violating LN by 2 units for Oi//3/3 decay is proportional to 
v^. This must be compared with the linear dependence in obtained when the x LN assignment is 2. 
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two to obtain the appropriate 6-fermion contact interaction 

= («7^(1 - 75)rf) (^7^(1 - 75)rf) e(l - 75)e'= , (24) 

where rrip denotes the proton mass and 

^3 = f^gl ■ 25 

This type of interactions has been aheady considered in the hterature [49], where hmits 
from the most sensitive experiments at that moment were derived. Since they have not been 
substantially improved, we will directly use the results in [49] from the Heidelberg-Moscow 
experiment corresponding to T1/2 > 1-9 x 10^^ years which yields 63 < 1.4 x 10^® at 
90% CL. On the other hand, experiments in the near future will be sensitive to lifetimes of 
the order of 6 x 10^'' years [14], i.e. a reduction factor on the coupling of roughly 0.05. Then, 
in order to Oz//3/3 decay be observable in the next round of experiments but still satisfy the 
present limits, we must require 

„ Next m^iii^vl n ^ . s 

8.75 X 10"^^ < r gee < 1.75 x 10"^ 90%C.L. , 26 

where rrip is the proton mass and the inequality with the superscript "Next" corresponds 
to the requirement that Ou/S/S decay will be observed in the next generation of experiments 
[14]. While the inequality without the superscript stands for the present experimental limit 
at the 90% C.L.. The conditions in Eq. (26) will prove rather restrictive because its range 
of variation is relatively narrow. In fact, reducing the lower limit will appreciable enlarge 
the allowed parameter region as discussed below. 

Thus, the above lower limit together with the requirement of perturbative unitarity (indi- 
cated by the "Pert" superscript in the inequalities below) or naturality, which bounds from 
above the product of couplings and VEVs /i^f^l^'eel, translate into an upper limit on the 
product of the scalar masses m^ia- These, however, are not precisely established because 
the perturbative bounds are in fact estimates that vary with the approach. In Fig. 2 we 
show the allowed — region, where m«; ~ and ^ in the limit of small 
mixing angle 9d, for = 2 (blue, darker) and 5 (orange, lighter) GeV (see Appendix A), 



There is a misprint in Ref. [49]. We thank the authors of this reference for providing us with the correct 
hniit on 63. 
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Figure 2: Projection on the — rriy plane of the allowed parameter space region, with 



and m-v 



On the left (right) we draw the region assuming perturbative unitarity (a LN 



breaking scale fi^ < 20 TeV). The blue, darker (orange, lighter) areas correspond to u^^- = 2 (5) GeV. 
The crosses stand for the reference point m^^^ = 10 TeV, m^j = 2 TeV (and = 2 GeV, /i^ = 
15 TeV, with g^e = 1 and g^^ = 0.001). 

respectively, assuming perturbative unitarity (left) and a maximum LN breaking scale 
(right). In the first case (see Eqs. (Bl) and (B2) in Appendix B) 



Pert , Pert 

\gee\ < V47r, /i« < 47r min(m«,^ 2) ; 



(27) 



whereas in the second one 



Pert 

< 



47r , fi^< 20 TeV . 



(28) 



It must be noticed that all (pseudo-)observables violating LN are proportional to fJ'n'V^- 
Hence, an increase in can be traded by the corresponding increase in /i^, and vice-versa. 
So, the orange areas in Fig. 2 can be also interpreted as the allowed regions for = 2 GeV 

Pert 

and < "ibiT min(mKi2) (left) and yU^ < 125 TeV (right). One may wonder at this point 
why we choose the bound of 20 TeV for /i^^ in Eq. (28); or equivalently, what is the effect of 
varying such a value. The answer is simple. The blue, darker region in the right panel of 
Fig. 2 disappears for /^^ ~ 8 TeV, which only reflects the narrowness of the range allowed 
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by Eq. (26), as required by our main working assumption that Oz//3/3 decay will be observed 
in the next round of experiments. The allowed regions in Fig. 2 are appreciably enlarged by 
reducing the lower limit in this equation. This can also be achieved by further increasing 

These areas are also bounded from below due to the non-observation of doubly-charged 
scalars; we can then assume 2 > 100 GeV, as discussed in Section VI. However, the 
enclosed areas in Fig. 2 are further reduced by a more stringent and subtle constraint. As 
we shall discuss below, bounds on LFV processes (see section IV) like — ^ e^fi^fi'^ banish 

to large values if the corresponding coupling product greO^fj, is sizeable, which is required 
because neutrino masses are proportional to Qab (see section V), and Qre and (7^^ must be 
large in order to accommodate the observed neutrino spectrum. Moreover, both and 

enter in the two-loop integrals generating neutrino masses, but these tend to zero in the 
limit m^/m^ — )■ 0. As a result, both scalar masses are constrained, but differently, by the 

—7- e^fi^fi" bound. The regions in Fig. 2 satisfy all experimental restrictions, including 
the upper bound in Eq. (26). The LHC will further reduce the allowed regions, mainly in 
the case of large LN breaking scale /i^. We also provide a "benchmark point", denoted by 
a cross in the figures, where all constraints are satisfied, and which we will use as reference 
throughout the paper. As can be deduced from this Figure the parameter space is rather 
constrained in this simple model when we require that the values of couplings and scalar 
masses stay natural, but one can think of other models within this class of theories where 
these constraints are significantly relaxed (at the price of complicating the spectrum through 
the introduction of additional scalars). 



IV. LEPTON FLAVOR VIOLATION CONSTRAINTS 

We will show in the next section that in order to obtain neutrino masses in agreement with 
experiment, the doubly-charged scalar Yukawa couplings Qab and the m^/m,^ ratio cannot 
be too small. In such a case some of the predicted LFV rates can become large enough 
to be at the verge of their present experimental bounds, especially for very rare processes 
like fi~ — )■ e~^e~e~ or /i~ — )■ e~7. Thus, we can use LFV processes to further constrain the 
model, and perhaps to confirm or exclude it in the near future. 

In this section we will briefiy discuss the most restrictive process — )■ e^e~e~, whose 
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tree-level amplitude is obtained by the single exchange of the doubly-charged scalar k. The 
corresponding branching ratio equals 

2 



BR(c ^ i+e 



c d . 



QabQcd 



BR(C ^ ^,-^^) , (29) 



2(1 + M 

where 6cd takes into account the fact that there may be two identical particles in the final 
state, as in our case, and 

1 _ cos^ 9d sin^ Op 

(If siv^Od <^ 1, the effective mass ~ since we, in practice, assume that m^^ 2 
are never very different.) Then, the current experimental limit on BR(yU^ — > e+e^e^) < 
1.0 X 10"^^ [1] translates into 

Ig^eglel < 2.3 X 10-^ (m,/TeV)^ (31) 

which is mainly a constraint on g^^ because g^f, must be relatively large if Oz//3/3 decay has 
to be observable at the next generation of experiments. 

Related processes provide weaker constraints. Thus, fi" — > e^7 proceeds at one loop and 
is suppressed by the corresponding loop factor, and similarly for /i — e conversion in nuclei. 
The bounds from /i+e~ < — > yU^e+ (muonium-antimuonium conversion) or muon-positron 
conversion, although tree-level processes, are also less restrictive (for a discussion of LFV 
processes mediated by doubly- charged scalar singlets in a similar model see [50, 51]). All 
these processes and the analogous ones involving r leptons, as well as the corresponding 
(anomalous) magnetic moments will be discussed in detail with more generality elsewhere. 

Here we shall be mainly interested in the interplay of a large Ou/SP decay rate and a real- 
istic pattern of Majorana masses, and for this purpose it is sufficient to show the restrictions 
on these (pseudo-)observables in simple SM extensions as the one at hand, and indicate 
which further processes may be within the reach of new experiments. In our case, the most 
restrictive process besides /i^ — )■ e^e^e" is — )■ e^fi^fi^, which will be discussed below 
taking into account the neutrino mass requirements. 



V. NEUTRINO MASS GENERATION AND 613 EXPECTATION 



In the model under consideration LN is not conserved when the couplings /i^, Ae , gab 
and Ye are non- vanishing. In this case there is no protection against the neutrinos acquiring 
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Figure 3: Two-loop diagram contributing to neutrino masses. 

Majorana masses m,^, which will then be proportional to all four couplings, are finite, and 
appear at the two-loop level, as explained in Section II and shown by explicit calculation in 
Appendix C. (If any of these couplings vanishes a conserved lepton number remains after 
spontaneous symmetry breaking and the neutrino masses will vanish.) These masses are 
generated by the non-renormalizable interaction aiiX'^L generated at two loops, when a 
and X develop VEVs; the corresponding coupling being mi,/v^v^. In contrast, the see-saw 
type II coupling iiX^L violates the Z2 symmetry and is forbidden to all orders. 

In Fig. 3 we draw one of the diagrams. Defining the neutrino mass matrix as usual, 

= -^TTZm^ul + h.c. , (32) 

we can write, taking <^ (see Eq. (A3) and Appendix A), 

('n..l.-^'nM, (33) 

where I^y is the sum of the (rescaled) loop integrals from the different graphs Jj, can be 
estimated in the mass insertion approximation with m^^-^ ^ mw- For instance, in this limit 
one of the contributions, Ji, corresponding to the diagram in Fig. 3 gives (neglecting the 
lepton masses in the denominator and assuming equal masses for the doubly and singly- 
charged triplet components) 

k^{k^ - ml)q%q^ - ml){{k - qf - ml) ' ^^^^ 



Again we note that neutrino masses, which violate LN by 2 units, are proportional to in our case. In 
contrast the dependence is linear in the estimate for the model in Refs. [21, 22], showing that x must be 
assigned LN 2, as in see-saw models of type II. Besides, the neutrino masses are in fact infinite in that 
particular point obscured because divergent diagrams were omitted. 
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In this approximation the full I^, is a dimensionless function oi y = {m^/irii^Y of order one, 
except for y — )■ 0, in which case it tends to zero as ylogy. In contrast with Ji which tends 
faster to zero, as (ylogy)"^, for vanishing y. is also bounded from above, going to a 
constant of order 2 for y — i- oo. A complete calculation taking into account the W-mass, 
as well as the corrections and the new scalar mass scales, is presented in Appendix C. 
A reasonable approximation is to neglect higher effects and take all triplet masses equal 
= iTLi_j = mA'i in fact, in the physical limit v^^ <^ v^^ k,2 is mainly the doubly-charged 
triplet component, u the singly-charged one, and A the imaginary part of the neutral triplet 
component. 

In our model the Yukawa couplings gab appear in the neutrino mass matrix, the Oi^/3/3 
decay amplitude and the amplitudes for the LFV processes, and this translates into rather 
stringent constraints on the allowed neutrino mass matrices (once one insists in dealing with 
a perturbative theory up to several tens of TeV). We now consider these constraints. 

Assuming /x^ ~ 10 TeV and < 2 GeV, and taking Ii,,\gee\ ~ 1? Eq- (33) gives 
I i''^'^)ee I ~ ^ 10^^ eV. How large can it be in general ? Using perturbativity limits 
(27) and 1,^ ~ 1, we get 



Pert ^ / mm T \ Next.Pert . , ^ 

< 1.6x10"^ f ] 1.6xlO-^eV, (35) 

where the upper limit is obtained by taking min(mKi 2) ~ 10 TeV (see Fig. 2, left). Alterna- 
tively, we can translate the limits on Oz//3/3 decay in Eq. (26) into bounds on (mj,)^^, but for 
large scalar masses this limit is less stringent than (35). In either case, | {rriu)^^ \ is typically 
less than 10~^. 

There is in addition a quite strong bound on [rriu)^^ from fi — )■ eee. Substituting Eq. (31) 
in the generic neutrino mass expression in Eq. (33) we find 

I im.) I < 2.3 X 10-^ ( -^^/. . (36) 

The constraint that a signal is seen in near-future Oz//3/3 decay experiments (left inequality 
in Eq. (26)) can then be used to eliminate \gee\'i in this way we obtain 

, , Next fh? vt Next,Pcrt ^ ^ ^ 

, I < 4.3/,^^^eV < 1.2xlO-5eV, 37 

where in the second inequality we used ly ~ 1, = 2 GeV, the naturality limit on /x^ 
(Eq. (27)), and smOu <^ 1. If we had used the LN breaking scale yU^ = 10 TeV and the 
production limit on > 0.1 TeV, then | {rriy)^^ | < 7 x 10^"^. 
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The final result of this phenomenological discussion is that both, | {m,y)^^ \ and | {friu)^^ \ , 
must be below ~ 10~"^eV, and this follows from requiring that (i) Oz//3/3 decay is at the 
reach of the next round of experiments, and (ii) that the theory is perturbative and free of 
unnatural fine tuning up to several tens of TeV. These limits could be somewhat relaxed: in 
the I (mjy)gg I case by making doubly- charged scalar masses larger, and in the | {m^)^^ \ case 
by allowing for a smaller \gee\- However, this is at the price of generating some tension with 
the naturality constraint in the former case, and spoiling the possibility of observing Oz//3/3 
decay induced by scalars in the near future in the latter one. There are additional but less 
severe bounds on the remaining Qab from other LFV processes; we discuss them below, when 
presenting the plots for the relevant neutrino mass (pseudo-)observables satisfying present 
experimental restrictions. 



A. Prediction for the third neutrino mixing angle 6*13 



The question now becomes whether it is possible to accommodate the observed spectrum 
of neutrino masses and mixing angles in this type of models once the above experimental 
constraints are imposed. In the following we will use the standard parameterization of the 
neutrino mass matrix [1, 48, 52, 53] in terms of 3 mass parameters, 3 mixing angles and 3 
phases: 

mv = UDyU^ , with = diag(mi, m2, ms) 



(38) 



and 



U 



C13C12 



C13S12 



-C23S12 - S23-Sl3Cl2e'" C23C12 - S23Sl3Si2e' 



S23C13 



v 



S23S12 - C23Sl3Ci2e' 



-S23C12 - C23Si3Si2e'" C23C13 



,401/2 



3*02/2 



v 



\ 



/ 



(39) 



where Sij = sin^y and = cos^y. A global fit to neutrino oscillation data gives (see, for 



instance, [54]) Am: 



(7.59+J;J°) X 10-5 eV^, Am 



10-3 eV', 



Q 019+0.017 2 



,+0.06 
-0.07' 



,2 
'13 



0.013 



31 



+0.007 



m. 



mr 



(2.50+^ 



+0.09 ^ 

16' 



X 



Neutrino oscillations are 



12 "-"""-O.OIS' "23 "•""_o.07' "13 """""" -0.005' 

not sensitive to the phases ai and 02, nor to a common mass scale which is conventionally 
chosen to be the lightest neutrino mass. 6, which appears multiplied by S13, is beyond present 
experimental sensitivity. The sign of Amg^^ is not presently known, and could be negative 



(known as inverted hierarchy), however, in this case 



m,. 



> 10 ^eV and cannot be 
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accommodated within our model; we will therefore consider only the normal hierarchy case 
Amg]^ > 0. Finally, recent data on electron neutrino appearance at T2K [55] and Double 
Chooz [56] experiments point out to a mixing angle ^13 different from zero. 

A possible way of identifying the allowed region in parameter space would be to first 
generate random values for masses, angles and phases within the 1 cr regions experimentally 
allowed in Eqs. (38) and (39), and obtain scatter plots for m,y. Then, using Eq. (33) we 
can solve for Qab, up to an overall factor oc fif^v'^I^/v^, and then find the values of fin,v^, 
mf^,m^ which respect the constraints discussed in the previous sections. The potential 
problem we face is due to the specific form of the neutrino mass matrix, which contains 
^aQab^b (see Eq. (33)), and is therefore suppressed for the first generations due to the light 
charged- lepton mass factors. To compensate this may require Qab to be too large to meet the 
bounds required by OuPP decay (Section III), LEV processes (Section IV) and perturbative 
unitarity (Appendix B). An alternative way to proceed is noticing that in practice (see (35) 
and (37)) we are asking if |(m,^)ee,e^| ~ is consistent with neutrino oscillation data (a 
question also of general interest not only within the model under consideration). These 
additional constraints will hold only within a limited region of the allowed neutrino masses 
and mixing parameters, which then implies that the type of models under consideration 
gives rather clear predictions about some of these parameters. 

In order to see how this comes about it is useful to go through a straightforward parameter 
counting exercise: mj, is a 3 x 3 complex and symmetric matrix specified by 12 real numbers: 
3 of these are unphysical and can be absorbed in re-phasing the neutrino fields, and 5 of the 
remaining 9 are measured (2 mass differences and 3 mixing angles, where we include 6*13). 
If we now impose {m^)ee,eti = 0, corresponding to 4 additional (real) constraints, only a set 
of points (or narrow regions, allowing for experimental accuracy) will be consistent. In fact, 
there may be no allowed values at all ! We have checked, however, that for each allowed 
choice of the experimentally measured parameters there is a unique solution for ai, ^2, S 
and mi satisfying all these restrictions. Eor example using the central values of the global 
fit given after Eq. (39) we find: 

ai = 0.65 , as = -2.32 , 5 = -0.78 , mi = 0.0036 eV, (40) 
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with 











0.59 + 0.58i 



m 



10 



,-2 







2.47 -0.2i 2.64 + 0.2i 



eV . 



(41) 



'U 



\ 



0.59 + 0.58Z 2.64 + 0.2i 2.12 - 0.21i 



This exercise can be repeated for different values of sin^ ^13, and the amazing resuh is that 
for the central values of Arrigi 31 and 23 thei'e are solutions only for 0.012 < sin^ 6'i3 < 
0.024, a range of values that roughly coincides with the result obtained by the global fits to 
present data and by recent T2K and Double Chooz experiments. To illustrate this result 
we present in Fig. 4 the sin^0i3 — 5 region allowed when \{rny)^e,efi\ = is imposed. The 
green, darker region is obtained when measured mixings and mass differences (except sin6'i3) 
are varied within 1 a, while the yellow, lighter one is obtained by varying them within 3 a. 
For comparison we also present the recent Double Chooz result [56] (sin^(2^i3) = 0.085 ± 
0.029 lb 0.042, adding statistical and systematic errors quadratically we obtain sin^ ^13 = 
0.022 ± 0.013) as a vertical band, while the cross stands for the reference point in Fig. 2. 

Of course, {my)^.^. and (m,^)e^ cannot be identically zero but small, ^ 10~^ eV. In fact, 
Qee must be different from zero and rather large in order to Oz//3/3 decay be observable (in 
this type of models {my)ee is small due to the huge suppression factor ml entering in its 
expression, not because gee is small itself). When {mv)ee and {m^)efx are allowed to vary 
within the model, with the other parameters staying within their 1 a range, sin^ 6'i3 is no 
longer bounded from above although the lower bound remains: 



In order to illustrate this behaviour we plot in Fig. 5 | {mu)e^\ as a function of sin^ 6*13, with all 
the neutrino masses and other mixing parameters varying arbitrarily within their 1 a limits. 
The red, darker region corresponds to |(mi,)ee| less than |(m,y)e^|. For comparison, we also 
plot the recent Double Chooz limit given above which is fully compatible with the T2K 90% 
C.L. interval 0.007 — 0.07 [55], and in agreement with current global fits (for instance, the 
one used in this paper [54] allows for 0.008 < sin^ ^13 < 0.020 at 1 a). Analogously, in Fig. 6 
we draw If function of the lightest neutrino mass mi. The red, darker region 

stands now for |(mj,)e^| < |(mj^)ee|- The previous figures are obtained from the neutrino 
mass restrictions only. But in our approach the model parameters are further constrained 
by the bounds from Oz//3/3 decay and LFV processes, which, as already emphasized, require 



sin^^i, > 0.008 . 



(42) 
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Figure 4: Allowed sin^ ^13 — 6 region for \{m,^)ee,efi\ = 0. The green, darker region is obtained 
when measured mixings and mass differences (except sin^is) are varied within 1 a; while the 
yellow, lighter one is obtained by varying them within 3 a. The middle dotted curve corresponds to 
the central values of the neutrino masses and mixings in the global fit performed in Ref. [54]. For 
comparison, we also draw the recent Double Chooz [56] 1 a limits, where statistical and systematic 
errors are added in quadrature (vertical dashed lines). The cross stands for the reference point in 
Fig. 2. 

\{^u)ee,e^i\ ~ 10^^ (see Eqs. (35) and (37), respectively). Then, 80^6*13 ^ 0.008 and 
0.002 eV ^ nil ~ 0.007 eV, as seen from Figs. 5 and 6, respectively. It then follows that a 
sufficiently precise measurement of sin ^13 can exclude the model being discussed, as stressed 
before. 

To conclude this section we derive the lower bounds, announced in Section III, on the 
scalar masses implied by the experimental limit on the — )■ e~^fi~fi~ branching ratio 
(< 1.7 X 10^®). In our model |(7refi'^^| niust be large in order to reproduce the observed 
pattern of neutrino masses. More precisely. 



\9eT9t^fi\ 



(2(27r)S^)^|(m.)^J|(m, 



lull Pert 

> 0.065 



TeV 



/,,min(mK, 2) 



(43) 



where the second inequality follows from Eq. (27) and the substitution of the other variables 
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Figure 5: \{m,y)efj,\ as a function of sin^^is. The green, lighter area which extends to all the panel 
is obtained varying the other measured masses and mixings within 1 a. The red, darker region 
satisfies the extra constraint |(miy)ee| < \irnu)eii\- The cross corresponds to the reference point in 
Fig. 2, and can be reached with either condition |(?Tij/)ee| > |(?TT'i/)e^|; whereas the darkest vertical 
1 a band stands for the recent Double Chooz result, summing errors in quadrature. We also draw 
the line |(m,y)e;x| = 10~^, which is the upper-limit estimate in this model. 
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Figure 6: |(m!/)ee| as a function of mi. Similarly to Fig. 5, but with the red, darker region instead 
satisfying the extra constraint |(miy)e^| < |(?7ij/)ee|- We also draw the line |(?7T.,y)ee| = 10~^, which 
is the upper-limit estimate in this model. 
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by their approximate values (see, for instance, Eq. (41) for the | (^i/)eT-/i^l estimates), in 
particular v^^ <2 GeV. Then, using Eq. (29) for — t- e^fi^fi^ and the experimental limit 
on its branching ratio one obtains l^'erfi'^t^tl < 0.007(mK/TeV)^, which combined with Eq. (43) 
yields > 1.2 TeV (see Fig. 2, left) We will present a detailed study of LEV processes 
in this type of models elsewhere. 

VI. COLLIDER SIGNALS 

Direct evidence for this type of models would be the discovery of the new scalars at a large 
collider. Doubly-charged scalars have fixed couplings to photons and are then produced at 
colliders with known cross sections. In addition their decay into leptons offers a very clean 
signal, which is particularly important at hadronic machines. Therefore, if doubly-charged 
scalars are light enough, they are very well suited for detection at colliders. 

In general, this type of scalars is assumed to be part of a weak triplet, and usually also 
acts as see-saw messenger of type II generating tree-level Majorana masses for the light 
neutrinos [46-48, 57] (see also [58]). Such triplets are then well-motivated on theoretical 
grounds, especially when considering LR symmetric models, and studies for searching the 
corresponding doubly-charged scalars at future colliders have been performed in the past [59- 
63] (see also [64, 65] for recent studies; model independent studies have also been carried out 
in the literature [66, 67]). The general conclusion is that the LHC discovery limit reaches 
masses over 600 GeV (for a center of mass energy of 14 TeV and an integrated luminosity 
of 30 fb^^) [64, 68]. However, the actual limits may be much better given the outstanding 
LHC performance, which almost matches the most favourable expectations for a CM energy 
of 7 TeV [69]. (See for a review [70].) 

Recently, first results from CMS have been presented at a CM energy of 7 TeV with an 
integrated luminosity of 0.89 fb~^ [71]. The analysis assumed a scalar triplet coupled to 
leptons, with 100% branching ratio to each leptonic channel. Nothing is seen, leading to a 
lower bound on the doubly-charged scalar mass of about 250 GeV if the main decay channel 
contain r leptons, and to about 300 GeV if they contain only electrons or muons. Weaker 
limits were obtained previously by LEP and the Tevatron. The absence at LEP of a pair- 

Note that we have used ^ 2 because is somewhat larger than to„ (see Appendix C). 
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production signal (e^e~ — )■ 7*,Z* — )■ kR) gives the constraint > 100 GeV [72-74] 
Limits on this type of scalars have been also derived using Tevatron data [76-78], leading 
to a limit > 100 — 150 GeV (depending on the details of the model). 

In our model the triplet does not directly couple to fermions, while the doubly-charged 
singlet does not couple to W pairs; however, triplet and singlet mix. Of the resulting mass 
eigenstates one, Ki, is mainly a singlet and decays dominantly to lepton pairs, while the 
other, K2, is mainly a triplet and has suppressed couplings to charged leptons. Both of them 
can be produced at LHC via the Drell-Yan mechanism (gg — )■ 7*, Z* — )■ ) with full 

strength. Since this is the main production process considered by CMS, the former limits 
apply directly to Ki for a non- negligible mixing: 300 GeV < . Limits on will be 
more difficult to obtain since the process qq — )• 7*,Z* — )• — )• W^W^W~W~ is much 
more complicated to deal with, due to its large backgrounds and the in general difficult 
reconstruction of several leptonic W decays. 

Notice that there are other production processes that are more model dependent. In par- 
ticular, the same interaction that induces Oz//3/3 decay and the decay of k into gauge bosons 
can mediate single-scalar production through WW fusion. The amplitude is proportional 
to the triplet VEV, f^, which is small, and to the singlet-triplet mixing, sin6'£), further sup- 
pressing this process. Nonetheless this could prove to be the dominant production channel 
at LHC [60] if f ^ > 1 GeV and 2 > 500 GeV. This is especially relevant for our model 
since both LHC and low-energy constraints require a relatively large v-^ as well as large 
scalar masses (see Fig. 2), all of which is in marked contrast to the constraints on triplet 
models with tree-level type II see-saw neutrino masses. A thorough study of the various 
possibilities is somewhat involved [79] and lies outside the scope of this paper; we will revisit 
these aspects of our model in a future publication. 

VII. CONCLUSIONS 

In this paper we have presented a simple model with a large Oz//3/3 decay rate into RH 
electrons through the exchange of the SM W boson and new heavy scalars x- this 

Single production via e^e^ — > nee, as well the u-channel contribution of k to Bhabha scattering have 
also been studied at LEP [74, 75], but the corresponding bounds depend on the unknown values of the 
Yukawa couplings. 
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model the light neutrino masses, Oz//3/3 decay and LFV processes have a common origin, 
which provides a simple description of these processes and also leads to a rather constrained 
parameter space. A Oz//3/3 decay final state with RH electrons also occurs in left-right (LR) 
models (see [80] for a detailed discussion), but generated by the exchange of new heavy 
neutrinos and gauge bosons Wr. In these models this process is a priori decoupled from 
the light neutrino mass generation, and hence the rate for Oz//3/3 decay can be large and 
at the same time the effective electron mass (mj,)ee small. Despite these differences both 
types of models have similar low-energy (below the electroweak scale) limits or, equivalently, 
they represent very different UV completions of two similar effective Lagrangians at scales 
below rnw- A general discussion of the different alternatives using an effective Lagrangian 
approach is presented in a companion paper [16]. 

The specific model we have discussed is one of a wide class of theories with similar 
phenomenology. For example, we choose to break LN explicitly by introducing one neutral 
scalar singlet (a), but, as discussed in Section II we could promote a to a complex singlet 
and insure the Lagrangian is LN invariant, and this symmetry spontaneously broken; while 
the neutrino phenomenology remains similar, it is simpler to discuss in our case. In general, 
this class of models require several new scalar multiplets to allow for non-vanishing scalar 
couplings to RH electrons, together with a chain of scalar couplings connecting them to the 
standard gauge boson triplet, and if desired, to the SM Higgs doublet in order to relate 
the new LN acting on the RH leptons to the ordinary LN associated with LH leptons. In 
such cases Oz//3/3 decay into RH electrons and the light neutrino Majorana masses are again 
related although at different loop order, as in our model. 

This model, which we considered in some detail, is based on an extended scalar sector 
containing three new multiplets: two isosinglets and one isotriplet, the physical spectrum 
contains two light scalars (with masses ~ "^^^o-) and 5 heavy scalars (with masses O(TeV)): 
2 doubly-charged, 1 singly-charged and 2 neutral. After writing the corresponding scalar 
potential and showing that there is a region of parameter space allowing for a local mini- 
mum along the desired direction, we have elaborated on their low-energy phenomenological 
implications. In deriving our quantitative predictions we required for the model to remain 
perturbative up to several tens of TeV. 

As repeatedly stressed, the new scalars mediate Oz//3/3 decay into RH electrons, LFV 
processes and, at two loops, generate Majorana neutrino masses. By requiring Oz//3/3 decay 
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to be large enough to be observable at the next round of experiments, we have derived 
lower bounds on the coupling to RH electrons gfge and upper bounds on the masses of the 
exchanged scalars. On the other hand, the stringent experimental limits on LFV transitions 
(particularly /x^ — )■ e^e^e^ and — t- e^yU^/x^) translate into upper bounds on the couplings 
Qab and lower bounds on the masses of the new particles involved, and these create some 
tension with the assumed Oz//3/3 decay rate. As a result, the model predicts that most LFV 
processes can be within the reach of the next generation of experiments, too. 

These same couplings and masses also enter the (two-loop) expression for the Majorana 
neutrino masses that have the very characteristic form, (mj^)af, oc rriaglfjmb, proportional 
to the scalar couplings to two RH leptons Qab and to the corresponding charged lepton 
masses rUa^b- Accommodating the observed neutrino mass spectrum and mixing parameters 
together with the other constraints is possible, but only for very restricted values of the 
lightest neutrino mass (not testable in neutrino oscillation experiments), and when the 
mixing angle 6*13 is constrained to lie within the preferred 1 a range from the global fit in 
Ref. [54], 0.008 < sin^^i3 < 0.020. This last prediction also goes far beyond our model, 
and holds whenever the (symmetric) neutrino mass matrix contains three very small entries, 
for example, rriee, men{= m^j_e) ~ 0, as in our case Some of these constraints can be 
alleviated by further extending the scalar sector, although at the price of requiring precise 
cancellations (that could be naturally enforced using further symmetries). 

The model considered presents a simple consistent extension of the SM exhibiting a large 
Oi//3/3 decay, but where the contributions to this decay generated by the neutrino Majorana 
masses are negligible. The neutrino masses themselves are predictable, in contrast with other 
proposals (see, for example [84]), and consistent with existing data. Note, however, that 
although our phenomenological approach and the constraints on the model mainly follow 
from requiring an observable Oi//3/3 in the next round of experiments, the non-observation 
of this decay and even a vanishing decay rate would be compatible with our analysis. That 
would be the case for g^e — > 0, value which would also ease the LFV restrictions without 
altering the neutrino mass predictions. As they are obtained assuming that (m,y)ee, which 
is proportional to gee-, too, is to a large extend negligible. 

A general discussion of the implications of texture zeroes in neutrino mass matrices can be found in 
[81, 82], and most recently in [83]. 
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Finally, we have also reviewed the collider limits on doubly-charged scalars. The excellent 
LHC performance should soon allow for the actual confrontation with the expected masses 
and couplings for the new scalars in the type of models studied here. We must, anyhow, 
be aware that once LHC settles the fate of the SM Higgs boson^^, the mass and couplings 
of the scalar doublet will be further constrained, implying further restrictions in the scalar 
potential which must be checked that can be satisfied. At any rate, the model studied here 
is one of a wider class sharing the main assumption in our analysis, that Oz//3/3 decay is large 
and decoupled from any mechanism providing tree-level neutrino masses, although the latter 
are generated through higher-order radiative corrections. 



Appendix A: Bound on the scalar triplet VEV 

As it is well known, the VEV of a triplet with hypercharge 1 gives a tree-level contribution 
to the p parameter, spoiling the successful SM (tree- level) prediction = 1. In general, 

- ^A^p? ' ^^^^ 

where the sum runs over the scalars of isospin Tj and hypercharge Yi with VEV fj. For the 
case under consideration, and assuming that the VEV of the triplet is much smaller than 
the one of the doublet. 

Then, the VEV of the scalar triplet contributes negatively to po; while the best p value 
obtained from a global fit to electroweak precision data is [1] p = 1.0008 J at la, 
and p = imOA+mil at 2a. Thus, vl/vl < 0.00035 at 2 a, implying < 3GeV for 

^ 174 GeV. Which is comparable to the bound derived from the global fit including 
explicitly the scalar triplet effects, < 2 GeV at the 90% C.L. [88]. 

However, a more complete analysis should also include the radiative corrections to the 

p parameter induced, for example, by the exchange of the scalar triplet, which can be 

positive [89]. For instance, in the triplet Majoron model [90] Ap = jJ-Gprn?, with 

27r^V 2 

In the very near future the LHC will be able to exclude a heavy Higgs with a relatively large significance 
[85, 86] but electroweak precision data do prefer a light Higgs (see, for instance, Ref. [87] and references 
there in), which shall require further for confirmation (or exclusion). 
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the mass of the doubly-charged scalar, cancelling partially the tree-level contribution. 
Since these contributions depend on the mass splitting of the triplet components, which in 
our model is not fixed, we will assume a conservative upper bound 



Appendix B: Constraints from naturality and perturbative unitarity 

The relevant parameter for decay and neutrino masses is the product of couplings 

and VEVs fi^v^Qee- It cannot be too large without leaving the perturbative regime and 
there are different arguments that can be used to set upper limits on its size (perturbative 
unitarity, naturality, etc.). 

Let us discuss the constraints from perturbative unitarity. Consider first the Yukawa 
couplings of doubly- charged scalars. Tree-level unitarity at high energy, s ^ ^kx2i i'^ 
ee — )■ ee collisions mediated by requires \gee\ < a/^tt- Similar bounds can be obtained 
from other channels. In order to be definite, we will demand 



Tree-level unitarity at high energy does not give useful information on dimensional param- 
eters like yU^ because amplitudes involving these couplings decrease with energy. However, 
it does not seem natural to have /if^ much larger than other dimensionful parameters in the 
model, like or m-^. In fact, one-loop diagrams involving the /i^ coupling give contribu- 
tions to or which are of order Sm^^^ ~ /i^/(47r)^. Therefore, it seems appropriate to 
require 



Limits (Bl) and (B2) also guarantee that the Hi decay width (to leptons and to K2, if this 
is light enough) is not too large as compared to its mass. 

One must be aware, however, that all these limits are estimates which depend on the 
naturality approach. Thus, although at the price of fine tuning, one might decide to fix the 
model parameters outside the range defined by these limits. Moreover, there can be model 
extensions where those values are natural. At any rate, we use Eqs. (Bl) and (B2) in the 
text to illustrate that the allowed regions in parameter space are at a large extent bounded 
if the perturbative theory must stay natural. 



< 5 GeV . 



(A3) 




(Bl) 




(B2) 
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Figure 7: Two-loop diagrams contributing to neutrino masses in the generalized gauge. 
Appendix C: Loop integrals for evaluating neutrino masses 

The evaluation of the complete two-loop contribution to neutrino masses, including m\Y 
effects in a general gauge, is a complicated task. Fortunately, we found a gauge in 
which the calculation simplifies enormously. Indeed, one can choose a gauge in which triplet 
and doublet charged scalars do not mix at all (the gauge-fixing Lagrangian is of the form 
Cgf = a\d ■ W'^ + bx^ + C(/)^p with the constants a, b, c chosen to cancel the aud the 

W^X~ y W~^(j)~ mixing terms). In this gauge the charged Goldstone boson becomes degen- 
erate with the physical charged scalar; and does not couple to fermions, whereas 
does not have derivative couplings to the doubly-charged scalars and gauge bosons. There- 
fore, the only diagrams contributing to neutrino masses are those depicted in Fig. 7. (The 
corresponding Feynman rules are obtained redefining the fields in Section II A adequately; 
in particular, in this gauge = and = 0^.) Following the notation in Eq. (33) we 
will have now I^y = Iw + I<i), where Iw and are the two contributions from the diagrams 
in Fig. 7. They read 

/„, = -2(irc fmi, cos* ffs J 4,2(^2 _ - ,ni-){{k - qf - - qf - mJJ ' 

( tA^ + A) + (1 - 4)' ,,2 tfl 2, ) ■ (CI) 



m 



w 



{Any /" ^'^ 



cos^ 61 J k'^{k'^ — mlj)q'^{q^ — m'^){{k — g)^ — m^J((A; — g)^ — m^^) ' 
where in we used the equality v^m\ = — Agf o-f ^ cos^ 6j (see Section II) to rewrite Ag in 
terms of itla- Note also that both doubly-charged scalar masses must enter symmetrically 
in the integrals to obtain a non-vanishing contribution. 

We have checked that in the decoupling limit {mw ^ '^ki,k2,i^,^ ^ Vcf,, implying 

171^2 = = cos ^5 = cos 6*/ = 1) we recover the results obtained in the mass 
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Figure 8: Contour plot for the loop integral ly as a function of the masses ~ and ~ rn^ 
for mw = 80 GeV and = ruA = ""t-kj- The cross denotes the reference point in the text. 

insertion approximation discussed in Section V. In general ly is a complicated function of 
m^y, mKi,K2.w,yi and the ratio v^/v^, which defines cos6'5 and cosOj. However, since <^ f^, 
we can safely neglect the corresponding corrections (in the limit of large m^^, the triplet 
VEV is small, and so is then the mixing of the triplet components with doublets and 
singlets; analogously, the contribution of doublet and singlet VEVs to the masses of the 
triplet components can be ignored). In this case 1^ becomes a function only of rn'f._^/'m'^ and 
m^^/m^ that we compute numerically. In Fig. (8) we plot the contours of constant 1^ as a 
function of the masses ~ and ~ ttL)^, where we have fixed mw = 80 GeV. From 
the Figure we see that is order one for a large region of the parameter space. Only when 
^ there is some suppression. The cross corresponds to the value of the reference 
point in Section III. 
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